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We consider the renormalisation of lattice QCD operators with one and two co- 
variant derivatives related to the first and second moments of generalised parton 
fN| ■ distributions and meson distribution amplitudes. Employing the clover fermion ac- 

■ tion we calculate their non-forward quark matrix elements in one-loop lattice per- 

<^ . turbation theory. For some representations of the hypercubic group commonly used 

in simulations we determine the sets of all possible mixing operators and compute 
\ the matrices of renormalisation factors in one-loop approximation. We describe how 

' tadpole improvement is applied to the results. 
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1> I. INTRODUCTION 

Many interesting observables in hadron physics, e.g. (moments of) generalised parton 
^ ■ distributions (GPDs) [l[ or distribution amplitudes, can be computed from matrix elements 
of local operators between hadron states. (For an extensive review of GPDs see Ref. 0], 
for distribution amplitudes see, e.g., Ref. [jjj].) GPDs, in particular, have attracted a lot of 
interest in recent years. They parametrise a large class of hadronic correlators, including 
e.g. form factors and the ordinary parton distribution functions. Thus GPDs provide a 
formal framework to connect information from various inclusive, semi-inclusive and exclusive 
reactions. Furthermore they give access to physical quantities which cannot be directly 
determined in experiments, like e.g. the orbital angular momentum of quarks and gluons 
in a nucleon (in a given scheme) and the spatial distribution of the energy or spin density 
of a fast moving hadron in the transverse plane. On the other hand, direct experimental 
information is limited and additional input is required to obtain a more complete knowledge 
of GPDs. One important source is lattice QCD, which can provide the relevant hadronic 
matrix elements @, H, @, 0, HI • 

Compared with moments of ordinary parton distributions, the specific difficulty in the 
treatment of moments of GPDs and distribution amplitudes lies in the fact that the required 
matrix elements of local operators are no longer forward matrix elements. In general, this 
circumstance complicates the pattern of mixing under renormalisation, and more extended 
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investigations become necessary. 

In a recent publication we have calculated the non-forward quark matrix elements 
needed for the renormalisation of quark-antiquark operators with two derivatives, which 
determine the second moments of GPDs, and we have discussed the mixing problem in 
detail. This calculation was performed in one-loop lattice perturbation theory for the Wilson 
fermion action. 

In the present paper we extend our perturbative calculations to improved fermions using 
the Sheikholeslami-Wohlert (clover) action [l(| for improving the vertices. We calculate 
renormalisation factors, but no improvement coefficients. While the general framework is 
of course the same as for Wilson fermions, the additional clover term in the action leads 
to a considerable complication of the calculations. A preliminary account of our work has 



already been given in Ref. [111 ]. Note, however, that a few misprints in Ref. [llj will be 
corrected here. 

Let us fix the notations used in our perturbative calculations. We work in Euclidean 
space and employ the Wilson gauge action together with clover fermions such that the total 
action is given by 

giatt __ 5< sw p _|_ Sw,g ■ (1) 



[x 



The fermionic part Ssw,f has the form [10 

t 
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CL 

— [^0)0 - ln)U x ,^(x + a/2) + ip{x + afi){r + 7/0^2,^ 0*0] (2) 



X,fl 

- i^f^E^^)^'^^° ve ^(x) 



written in terms of dimensionful massless fermion fields ip(x). Here a denotes the lattice 
spacing and the sums run over all lattice sites x and directions /i, v. All other indices are 
suppressed. F^ l ° VCT is the standard "clover- leaf" form of the lattice field strength (see, e.g., 



Appendix D in Ref. [12|). The coupling strength of the improvement term is given by c sw . 
The link matrices U X)fJi are related to the gauge field A^(x) by 

U x>fl = exp [igaA^x)] , A^x) = T c A^(x) , (3) 

where g is the bare gauge coupling and the T c are the generators of the SU (3) algebra. The 
gauge action for the gluon field A^x) is 



X,[l<V 



^Tr (U x ,, v + C/t J 



(4) 



with 

U x ,fiv — U x ^U x j ra f L ^U x+aj> fl Ul l/ . (5) 

In a perturbative calculation the investigated operators are considered between off-shell 
quark states. Our calculations are performed in Feynman gauge, the final numbers will be 
presented for the Wilson parameter r = 1, but for arbitrary values of c sw . 

On the lattice the operators are classified according to the irreducible representations 
rjp of the hypercubic group H(4) (for the notation see, e.g., Ref. [13]). Here / denotes 
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the dimension of the representation and k labels inequivalent representations of the same 
dimension. In addition our operators will be chosen such that they have definite charge 
conjugation parity C. 

Using clover fermions, quark-antiquark operators with one covariant derivative have been 



discussed in |14j for forward matrix elements. Because of the constraints imposed by charge 
conjugation invariance no additional mixing has to be considered in non-forward matrix 
elements and the renormalisation constants given there can be taken over to the case at 
hand. We will present them again for completeness and add the corresponding results for 
operators involving [7^,7^], which were not considered in Ref. HI 



In j9| the renormalisation procedure for the case of non-vanishing momentum trans- 
fer has been discussed in detail. We will not repeat this discussion here. The matrix of 
renormalisation and mixing coefficients Zij(a/j) relating regularised lattice vertex functions 
Vf(j/,p, a, <7r) and MS renormalised vertex functions Tf{p' ,p, fi, g^.) is defined such that 





N 



r?(p',P,fi,gji) = Z^ia^^Zijia^Tfip 1 ',p,a,g R ) (6) 

i=i 

with the quark wave function renormalisation constant Z^. Here p (p f ) denotes the mo- 
mentum of the incoming (outgoing) quark, the renormalisation scale is /1, the renormalised 
coupling is denoted by g^, and N is the number of operators which mix in the one-loop 
approximation. 



II. OPERATORS AND MIXING 

We consider operators with up to two covariant symmetric lattice derivatives D = D — D 
and external ordinary derivatives d. In the operator symbols, the derivatives are indicated by 
superscripts D and d. Note that matrix elements of operators which are ordinary derivatives 
of other operators are simply given by the matrix elements of these other operators multiplied 
by the appropriate product of components of the momentum transfer. 

The quark-antiquark operators with one derivative are given by 

^h„5„v, (7) 

i *~* 

-tp^j 5 D^ , (8) 

7^[7m,7^]A>, (9) 

^■,Ofc,7*W0 • (10) 

Operators such as ([9]) involving [7^, 7„] are of interest for tensor GPDs as well as for transver- 
sity and we call them transversity operators. They are antisymmetric in the indices \x and 
v. For non-chiral fermions, operators ([§]) and ffTUl) contribute as lower-dimensional operators 
to mixing in certain operators which determine second moments of GPDs. 



O 



D 



o 



5,D 



o 



T,D 



T,d 



i 



As operators with two derivatives we consider 

0°Z = --3l,D u DW 
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0%L = -f„[^D^) , (11) 
OIL = ~\d»d u (^) 



and 



0% D f = -- A ^1,D U D^ 



1 



0%l D = --.dA^^DM , (12) 



4 

o%2 = --Ad. (^7,75^) 



4 

We include here also transversity operators with two derivatives: 



O 1 ^ = —^b„ lAD.D^ , O^Z = (ife, 7«W ■ (13) 

A detailed description of the operators with two derivatives, their representations and poten- 
tial mixings is given in To define the various operators we use the following short-hand 
notations: 

0... {uiU2 } = - {0... UlV2 + 0... V2Vl ) , (14) 

0\\viV2V3\\ Ou\V2V3 Ov\V3U2 Ov3VlV2 ^U3U2V± 2 (D U2 U3Vi 2 Ou2VlU3 ) (1^) 

For the first moments we choose the following representations and operators (for a more 
detailed discussion of the transformation under H(4) see fl3|): 

Operator Representation C 

Ofu } +1 

og - 1 (cfi + og + o 3 D 3 ) n (3) +1 

{ 5 i? } "I (18) 



rf -1 



Off - I + ^ + Olf) 

T,D (8) 

Q T,D _ T,D (8) x 

^<{122)> ^{(133» '1 ^ i 



All operators in (fl8"|) are multiplicatively renormalisable. These representations exhaust 



all possibilities for the twist-2 sector in the continuum. Note that in Ref. [ll| we have 



erroneously assigned t\ and C = — 1 to an operator belonging to t^ ' , C — +1 



(8) 
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Let us now turn to the second moments and the corresponding twist-2 operators. In the 
unpolarised case we consider the following sets of mixing operators: 

T?\C=-l 



-1 

o 2 = 

/n _ ro- - x (n DD 4- n DD \ 

U 3 - <^«114» - 2 l U ({224)> + U ((334»J 

C 4 = Cf(ii4)) - - (Cf{224)> + Cf(334») > ( 2 0) 

e> 6 = e> 
c 7 = i 

Os = _ 

There is one more representation, t[ 4 \C = — 1, giving twist-2 operators. However, even 
in forward matrix elements the corresponding operators mix with operators whose dimension 
is smaller by two. Therefore they are rather unsuitable for numerical simulations and will 
not be discussed any further. 



\DD 


1 (n DD 4. 


n DD \ 


{114} " 


1 in" 4- 

~ 2 l U {224} + 


^{334}^ > 


\DD 
«114» 


1 (o DD 

- 2 l U «224)> 


+ L/ «334» 


,9.9 

«114» 


1 (O 99 

2 l U «224)> 


+ U «334» 


,5,(91) 






1 1 213 1 1 ' 












«213)> 






,5,-D-D 






1 1 213 1 1 ' 






{T,d 
411 


\ (on + o 


T,d\ 
433 J • 



In the polarised case we consider 



^{124} ' ^{124} ' \ Zl ) 



A 8 \c 



n 5 — n 5 ' DD — ( n 5 ' DD _i_ n 5 ' DD \ 

U l ~ U {114} 2 \ {224} + U {334} J ' 

sr\b srS,dd 1 / s/yb,dd _|_ sr\5,dd \ 

U 2 — U {114] 2 V i 224 } ^{334} J ' 

— n 5 ' DD ^ ( n 5 > DD _l n 5 > DD \ 

U 3 ~ U ((114)) 2 V «224» + ^«334» ) > 
4 _ ^<{H4)> 2 V«224» (<334» ^ ' 

/n5 _ /-oar* 

^5 — ^1121311 ' 

^6 — ^((213)) > 

U 7 — ^1121311 > 

/rj5 ey^^ o^o^'^ cy^'^* 

^8 ~~ ^123 ZLy 231 ^132 ■ 
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Again, there is one more representation giving twist- 2 operators. Here it is r| , C — +1, 
but the corresponding operators suffer from similar mixing problems as the operators with 
t\ , C = — 1 in the unpolarised case and will be omitted in the following. 

Finally the transversity operators: 



QT _ Q T,DD qT _ T,dd f2 o\ 
^1 ~~ ^411231 ' ^2 — Ly 4U23> ' \ ZO ) 



T<?\C 



n r _ _n T > DD _l n T ' DD — n T ' DD _l n T ' DD _ nn T ' DD 

U 3 — °1{133} °1{144} °2{233} °2{244} ZU 3{344} ' 

n T — -n T ' 99 j- n T < 99 _ n T ' 9a _l r) 5 "' 99 _ on T >" <oa\ 

^4 — Ly i{144} L/ 2{233} L/ 2{244} ZLV 3{344} ' 



n T — n T ' DD _l n T ' DD _ n T ' DD _ n T < DD 

U 5 — L/ 13{32} ^S^l} L/ 14{42} L/ 24{41} ' 

n T - r> T > 9,9 _l r> T > 9,9 _ ro T ' ,9a _ n T ' 99 (ok\ 

^6 — ^13{32} ^23{31} ^14{42} ^24{41} ■ l Zd J 

The representations fl23|) . (1241) and fl25l) exhaust all possibilities for transversity operators 
of twist 2 which have a mixing matrix of size 2x2 only. Other representations have more 
complicated mixing patterns. 



III. ONE-LOOP CALCULATION 

We calculate the matrix elements of the operators in one-loop lattice perturbation the- 
ory in the infinite volume limit following Kawai et al. [15|. Details of the computational 
procedure, in particular the Feynman rules, are given in [9|. 

A. First Moment 

Since mixing is absent we need no matrix indices for the renormalisation constants and 
use the general form 

ZM = 1 - ( 7 ln(aV ) + B(c sw )) , (26) 

where the finite piece B depends on c sw and the (one-loop) anomalous dimension 7 is given 
by 

/ 8/3 for rf ,T<ViVf 
H 3 for rf\r? ' ™ 
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For the operators ()18p we have (results for r 3 , r| , r 4 and r 4 are taken from 14]) 
Representation B(c sw ) 



r (6) 


1.27959 - 


3.87297 c sw 


- 0.67826 c 2 sw 


(3) 
r l 


2.56184- 


3.96980 c sw 


- 1.03973 cL 


r (6) 
'4 


0.34512 - 


1 .35931 


- 1. 89255 cL 


r (3) 
'4 


0.16738 - 


1.24953 c sw 


- 1.99804 cL 


r (8) 
'2 


1.25245 - 


3.10180 c sw 


- 1.59023 


(8) 


0.52246 - 


2.99849 c SU) 


- 1.46224 cL 



(28) 



B. Second Moment 

We write the matrix of renormalisation constants in the generic form 

Z\?M = % - M«V) + 4 m) (c„)) , (29) 

with 1 

dW/„ \ _ p(0,m) , R (l,m) R (2,m) 2 /orA 

The superscript (m) with m = 1,11 distinguishes the realisations I and II of the covariant 
derivatives, which are explained in Appendix A of Q. In the first case the momentum 
transfer "acts" at the position x associated with the operator, where we define for one 
covariant derivative 

0{x) = (tfD^ (x) = —(^(x)U Xtfl ijj(x + ap) - *P(x)Ul^ a ^{x - ajj) 

+ $(x - afi)U x _ afl ^(x) - $(x + ap)Ul 4l ^{xYj . (31) 

We have set the Dirac matrix in the operator equal to the unit matrix for simplicity. Reali- 
sation I leads to 

(V^B^) ( \) = Y a E + - + aft)Ul^(x)] + e^+^j .(32) 

X 

Alternatively (realisation II), q can be applied at the point half way between the ^ and ijj 
fields: 

/_<-» \ ( n ) 1 

W (?) = - E WO^O* + «/0 - ^ + a/i)^(x)] . (33) 



1 B'lj''"^ has been denoted by c-"^ in 



8 



In the case of two covariant derivatives we have for realisation II 

/_<-<- \ ( /7 ) 1 /- 

( V'-D/i-Di/V'J (q) = — } j {^{x)U XilJ U x+a(i ^{x + ap + oi>) 

-^0 + ai>)U x+ai>tlM Ul +a(ltU i/;(x + ap) 
—ip(x + ap)UljJ x ^{x + az>) 

and realisation I is obtained from (1341) as 

^3^) I] (?) = cos (^) cos (^) (^B M B^) (//) (g) • 
We get the following results: 



(34) 



(35) 



In this case we have the mixing operators ( JT9l) . The corresponding 2 x 2-mixing matrices 



arc 



7 



5(0,/,//) 
5(2,/,//) 



25 _5 
6 6 



11.56318 0.02414 \ 
20.61780 J 

2.89800 0.25529 
-4.74556 

0.98387 -0.01557 \ 
-0.54317 ) 



(36) 
(37) 
(38) 
(39) 



The matrix B^ 1 ' 11 ^ shows a rather small coefficient for the mixing between the operators 

n DD anrl O" 
U {124} allCL U {12 

where /1 = 1/a. 



Of^ and Cff 24} . Thus it may be justified to neglect the mixing in practical applications, 



Q 1 [rf\C = -\) 

The operators mixing with 0\ are given in fTSOj) . First we consider the operators of the 
same dimension Oi, . . . ,Oj. To one-loop accuracy the operator O7 does not contribute, 
because its Born term vanishes, and we have to consider the following mixing set: 

{O l7 2 ,0 3 ,0 4 ,0 5 ,0 6 }. (40) 

The anomalous dimension matrix is 

/ f -§ o\ 




7 



I -| 




2-2 

\ 00 o-l 1/ 



(41) 
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and the finite parts of the mixing matrix are given by (in doublets the upper number belongs 
to type I, the lower to type II of the realisation of the lattice covariant derivative) 



/ 



■12.12740 


3.30605 







1.49127 \ 
-2.73669 J 
20.61780 

-8.01456 \ 
18.18411 / 


3.26440 
3.26440 



0.36848 


-14.85157 







-0.41595 \ 
0.99336 ) 




0.01562 0.14983 




\ 



-0.92850 0.73802 



4.43061 

-4.30228 
20.61780 

0.35008 0.01491 

0.00497 0.36003 / 



(42) 





f 0.21269 \ 


0.03276 


2.92169 




^ 0.68643 J 







-4.74556 







( 0.76570 \ 




0.33335 


-2.15228 




V 0.05510 J 
















1.44106 








1.44106 






-0.01492 
-0.17283 







-1.20652 
-0.96966 
-4.74556 




0.01878 -0.05696 



1.75814 -2.29837 











■1.64790 -0.86576 



-0.28859 -2.22507/ 



(43) 



/ 



_0(2,/,//) 



V 



-0.98166 


0.37050 







-0.07815 
-0.10117 
-0.54317 

-0.55068 \ 
0.21545 J 


1.41570 
1.41570 



-0.02914 


-1.70741 







0.03475 
0.04243 







0.37132 \ 
0.11594 J 
-0.54317 




-0.00999 0.00688 






-0.44295 


-1.70334 
0.18921 



0.10328 


0.56763 
-1.32493 / 



(44) 

The matrices B^' 1 ' 11 ' show sizeable coefficients for the mixing of 0\ with other operators, 
especially with O2 containing two external ordinary derivatives. 

There is also a possible mixing between 0\ and the lower- dimensional operator (9 8 in 
( 12"U|) . Indeed, we find in the one- loop approximation that the vertex function of 0\ contains 
a term oc 1/a: 



o 1 



-part 



9 2 rC f 
16vr 2 



-0.51771 - 0.08325 c,„, - 0.00983 d 



1 



O 



Born 



(45) 



10 



This mixing leads to a contribution which diverges like the inverse lattice spacing in the 
continuum limit. Thus the perturbative calculation of the mixing coefficient is not reliable 
and the operator O s has to be subtracted non-perturbatively from the operator 0\. 



q5,DD ( J4) 



{124} 



4 V ,C = +1) 



In this case we have to consider the operators in (I2ip . The anomalous dimension matrix 
is given by (1361) . the finite contributions are collected in the matrices 



5(0,/,//) = 



-12.11715 0.16673 \ 
15.79628 / 

-1.51925 0.00505 \ 
0.24783 J ' 

-1.71846 0.00711 \ 
-2.25137 / 



(46) 

(47) 
(48) 



As in the case of the operator 0?^, above, the mixing coefficient is rather small. 
0\ (t 2 {8) ,C = +1) 



First we discuss the mixing of operators of the same dimension in ( 1221) . The set of 
contributing operators is found to be 



{olololololot} . 



(49) 



As in the case of the operator 0%, one operator - here Oj - does not contribute to mixing 
in one-loop order. The finite contributions are 



/ 



-12.86094 


3.42196 







-2.06532 
1.48943 
15.79628 

15.82073 
-7.30020 



-8.91237 
-8.91237 



0.34900 


15.35920 







0.85381 
-0.33110 



0.05113 0.05942 



\ 







-5.16392 \ 
2.54306 ) 
15.79628 



0.17014 -0.94314 











0.95969 -0.95969 
-0.31990 0.31990 / 



(50) 
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V 



-1.49316 


0.09099 





-1.68673 


0.16581 







0.03027 \ 
-0.17480 J 
0.24783 

-0.99704 
1.07971 



13.26724 
13.26724 

0.03501 \ 
0.19704 ) 
-2.25137 



0.28263 
0.13256 



-1.41570 
-1.41570 



0.00750 


-2.30129 







-0.00612 


-1.36546 







-0.03022 \ 
0.09858 J 


1.25511 
0.56286 
0.24783 



-0.01275 \ 
-0.06677 J 


-0.30614 \ 
-0.25600 J 
-2.25137 



-0.01290 0.03622 



-0.13420 -0.48232 











-3.27954 1.78029 
0.59343 -2.09268 / 



0.01000 -0.00692 



\ 



(51) 



0.44299 -0.10322 











-0.97445 -0.88857 
-0.29619 -1.56683 / 



(52) 

The anomalous dimension matrix is the same as for the operators without 75, see ( l4"Ti) . 
Again, some of the mixing coefficients are non- negligible. The mixing between 0\ and Of 
is present also in the forward case. 

We also find mixing with a lower-dimensional operator, in this case it is the operator Of 
in ( |22|) . The corresponding contribution in the vertex function of 0\ reads 



-part 



skc t 

16tt 2 



(0.25231 - 0.02507 c sw + 0.01046 c 2 sw ) - O 



5, Born 



(53) 



Of(r 2 (3) ,C = -l) 
For the operators ( 1231) the matrix of anomalous dimensions is given by 



7 



and the finite contributions are 



13 _2 

3 3 

1 



-11.54826 0.21894 \ 
17.01808 

-2.41077 -0.05383 
-3.91333 

-1.51175 -0.00614 \ 
-1.97230 J 



(54) 

(55) 

(56) 
(57) 
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The operators (124)) have the same anomalous dimension matrix (1541) as the previous case. 
For the finite pieces we obtain 

B {o,i,ir) = ( -H-86877 0.27533 \ 
" V 17.01808 J ' 



B (i,i,ii) = ( "2.30651 -0.01831 \ 

V -3.91333 J ' 1 J 

s (V,/i) = ( -1.34908 0.01726 \ 

V -1.97230 J ' V ; 



or(r 2 (6) ,c = -i; 



We find for the case (|2"5"|) the finite mixing contributions 

(0.J.JJ) = f -11.74773 0.23797 \ 

" V 17.01808 J ' 1 J 

(WJ) = / -2.36201 -0.04490 \ 

V -3.91333 / ' 1 ; 

B ( W i) = / -1-45084 0.00898 \ 

V -1.97230 J { ' 

with the anomalous dimension matrix (15 



IV. TADPOLE IMPROVEMENT 

It is well known that many results of (naive) lattice perturbation theory do not agree 
very well with their counterparts determined from Monte Carlo calculations. One main 
reason for these discrepancies is the appearance of gluon tadpoles, which are typical lattice 
artefacts. They turn the bare coupling g into a poor expansion parameter. As a remedy the 
so-called tadpole (or mean field) improvement has been proposed [l6|, a rearrangement of 
the perturbative series making use of the variable Uq, the fourth root of the measured value 
of the plaquette, 

u ={-±-TrUn)*. (64) 

Its value depends on the coupling g 2 = 6//3 where it has been measured. 

In case of mixing the tadpole improvement works as follows. Scaling the link variables 
with m 

U,{x)=u (EM\^ UQ u^ x ) ( 65 ) 
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u = l-^ n 2 + 0(g*). (67) 



one finds an expression for the vertex function Tj of an operator Oj containing rij covariant 
derivatives which is of the form 

T 3 (U,(x)) = u^T 3 (U,(x)). (66) 

Tj(Ufj,(x)) = u T^U^x)) is expected to have a better converging perturbative expansion, 
which is obtained by inserting the expansions of ii " 3 and Yj{U^{x)). To one- loop accuracy, 
uq is given by 

16 7T 2 

Note that at the one-loop level we do not have to distinguish between g and #r in the 
perturbative expressions. The exponent of Uq depends on j because in general the mixing 
operators have different numbers of covariant derivatives (see, e.g., ( I19p or (I20p ). An external 
ordinary derivative (d) does not provide a factor of u . Taking into account the mean field 
value for the wave function renormalisation constant for massless Wilson and clover fermions 

Z? F = u (68) 
the mean field Z factor for each operator Oj contributing to the mixing reads 

Zf F = u~ n > , (69) 
and the tadpole improved matrix of renormalisation constants is given by 

rvMF / „2r< 



(70) 



Additionally, we replace the parameters g and c sw by their "boosted" counterparts 

9ti = 9 2 «o 4 , c T s l = c sw u 3 . (71) 
Combining (129~|) . (1701) and (1711) we obtain for the tadpole improved matrix in one- loop order 

l ~ ni ( ~ 4# f 7« V) + <' (m) (<£)) 1 • (72) 



Let us exemplify the impact of tadpole improvement in some typical cases. We choose 
\i = 1/a, c sw = 1 + 0(g 2 ) and u = 0.8778 corresponding to quenched calculations at j3 — 6. 
For operators with one covariant derivative the tadpole improvement procedure is rather 
simple. Because n\ = 1 the only effect consists in replacing in (1261) and (125]) c sw by and 
g by gxi- For the operators belonging to the representation r!f\ e.g., we get 

Z = 1.028 -> Z T/ = 1.023. (73) 

The operators for the second moments of GPDs are a bit more involved. First we consider 
the simple mixing 0^ l2A ^ <-> C^ 24 , (J2TJ) . Without tadpole improvement we obtain the mixing 
matrix of renormalisation constants as 



1.12965 -0.00151 
0.88354 



(74) 
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The tadpole improved result is 



7 Ti 



1.20501 




-0.00216 
0.81458 



(75) 



It is instructive to compare the one-loop corrections for the renormalisation constants, 
B^" (c 8W ) for the unimproved case ( |29i) and (c^) for the tadpole improved case 

2]). With the parameters given above we get for the operators ( |2~TT) 



B 



-15.35486 0.17889 \ 



and 



B 



TI 



13.79274 

-4.06129 0.17340 
5.06434 



J 



(76) 



(77) 



Equations (IT6]) and f!77|) show that the one-loop corrections on the diagonal have been 
reduced significantly. This is in accordance with the aims of tadpole improvement. 

The same procedure can be applied to the more complicated set Ox, ■ ■ ■ , O e (t[ 8 \ C = — 1) 
from (T20]) . We obtain for the unimproved case the mixing matrix of renormalisation constants 
(lattice covariant derivative type m = I) as 

/ 1.13535 -0.013727 -0.00314 0.00334 -0.00021 -0.00084 \ 



Z (i) 



\ 





-0.02823 






0.87057 
0.06585 


-0.05168 
-0.05168 





1.15799 














-0.03036 -0.00327 



0.87057 







1.02534 
0.00080 





0.01230 


0.00239 
1.02693 J 



(78) 



The tadpole improved result reads 



Z TI,(I) 



( 1.21508 -0.01997 -0.00611 0.00512 -0.00034 -0.00163 \ 

0.78680 

-0.05265 0.09673 1.25616 -0.04724 -0.00082 0.01094 

0.78680 

-0.06100 1.02195 0.00442 

V -0.06100 0.00147 1.02490 / 



(79) 



The renormalisation and mixing matrices for all other cases can be treated analogously. 
They show a similar behaviour: The 2x2 matrices exhibit rather small mixing coefficients, 
the other 6x6 problem fT22|) is almost identical to the example discussed above. 



V. SUMMARY 

In this paper we have considered quark-antiquark operators needed for the computation 
of the first two moments of GPDs and meson distribution amplitudes within the framework 
of lattice QCD. In one-loop lattice perturbation theory we have calculated the non-forward 
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quark matrix elements of these operators employing clover improved Wilson fermions and 
Wilson's plaquette action for the gauge fields. From the results we have evaluated the 
matrices of renormalisation and mixing coefficients in the MS-scheme. 

For the operators with only one derivative (relevant for the first moments) we could take 
over the numbers obtained for the first moments of ordinary parton distributions. The 
results for the second moments generalise the numbers calculated with Wilson fermions 
The general conclusions concerning the mixing properties remain unchanged. 

If there is only mixing between one operator with two covariant derivatives D and one 
operator with two external derivatives d the mixing coefficient turns out to be rather small. 
In the two cases (120]) and ff22l) with eight potentially mixing operators the mixing effects are 
more severe. Moreover, taking 0\ from (1201) or 0\ from (122]) as the operator to be measured 
in a numerical simulation we find in each case mixing with a lower-dimensional operator 
(cf. f|4~5"j) and fl5"B~l) ). This could lead to difficulties, because 1/a effects are hard to get under 
control. For overlap fermions, however, these mixings with lower-dimensional operators of 
different chirality would be absent. 

Additionally, we have discussed tadpole improvement with special attention to mixing 
operators. We have given the general prescription for mean field improvement and have 
shown how it works for selected examples. 
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